aspect = 4.0/2.2
scale = 5


#set output 'ss.gif'





set key top right
set grid lt 1 lc rgb "#EDEDED"


set xlabel "x [sigma]"
set ylabel "y [sigma]"

dr =  2**(1.0/6)
set xtics 1
set ytics 1

_log(a, b) = log(b)/log(a)
a(z) = 1.0/(z+1)
	
n = 1300.0
si = 400.0
z0 = 127.0
zi = 1.0
p = _log( si/n,  (a(zi)-a(z0))/(1.0-a(z0))  )

set terminal gif animate delay 5 font "Arial" 12 enhanced size 650, 800
set output 'ss.3d.gif'
i = 0
unset border
unset xtics
unset ytics
unset ztics
unset grid
unset xlabel
unset ylabel
unset zlabel
while (i<=20) {
	i=i+5
	#set title sprintf("p.size() = 4096, step = %i", i)
	
	splot [-9:9][-9:9][-9:9] sprintf('ss/%03i.txt', i)  u ($1/dr):($2/dr):($3/dr) with p lc 1 ps 1.9 lt 7 t "";
	
	print "i=", i
	
}

set terminal gif animate delay 5 font "Arial" 12 enhanced size 600, 600
set size square
exit
set output 'ss.gif'
i = 10
while (i<1990) {
	set title sprintf("p.size() = 4096, step = %i", i)
	plot [-9:9][-9:9] sprintf('ss/%03i.txt', i)  u ($1/dr):($2/dr) w p pt 7 t ""
	
	
	print "i=", i
	i=i+10

}
